The impact of correlated projections on weak lensing cluster counts 



o 
o 

(N 
o 

Q 



O 
U 

o 



O 

o 



X 



Laura Marian^, Robert E. Smith^, Gary M. Bernstein"^ 

^ Argelander-Institut fiir Astronomie, Universitdt Bonn, Bonn, D-53121, Germany 
^Institute for Theoretical Physics, University of Ziirich, Zurich, CH 8037, Switzerland 
'^Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA 

Large-scale structure projections are an obstacle in converting the shear signal of clusters de- 
tected in weak-lensing maps into virial masses. However, this step is not necessary for constraining 
cosmology with the shear-peak abundance, if we are able to predict its amplitude. We gener- 
ate a large ensemble of A'^-body simulations spanning four cosmological models, with total volume 
Kot ~ Gpc]"^ per model. Variations to the matter density parameter and amplitude of fluc- 

tuations are considered. We measure the abundance of peaks in the mass density projected in 
sslOO Mpc slabs to determine the impact of structures spatially correlated with the simulation 
clusters, identified by the 3D friends-of-friends algorithm. The halo model shows that the choice 
of the smoothing filter for the density field is important in reducing the contribution of correlated 
projections to individual halo masses. Such contributions are less than 2% in the case of the opti- 
mal, compensated filter used throughout this analysis. We measure the change in the mass of peaks 
when projected in slabs of various thicknesses. Peaks in slabs of 26 Mpc and 102 Mpc suffer 
an average mass change of less than 2% compared to their mass in slabs of 51 Mpc. We then 
explore the cosmology dependence of the projected-peak mass function, and find that, for a wide 
range of slab thicknesses (< 500 Mpc), it scales with cosmology in exactly the same way as the 
3D friends-of-friends mass function and the Sheth-Tormen formula. This extends the earlier result 
of Marian et al. [S^]. Finally, we show that for all cosmological models considered, the low and 
intermediate mass bins of the peak abundance can be described using a modified Sheth-Tormen 
functional form to within 10%-20% accuracy. 



I. INTRODUCTION 

For more than a decade, weak gravitational lensing 
(WL) has been considered a powerful probe for testing 
the cosmological framework. The basic idea of gravita- 
tional lensing is that on its way to the observer, the light 
from distant sources is deflected by the intervening mat- 
ter along the line of sight. On average these deflections 
are very small, so the result in most cases is that a circu- 
lar source galaxy looks slightly elliptical, i.e. the image 
of the source is sheared. Images can also be magnified. 
For a quantitative idea, the shear is typically of about 
2%, hence the name of WL. By measuring the shear, we 
can map the 3D matter distribution of the universe in an 
unbiased way, independent of baryonic matter tracers. 
Indeed, the largest surveys to date aimed at measuring 
cosmic shear, from CTIO (Jarvis et al. llaJlTl ) and more 
recently the CFHTLS (Hoekstra et al. [lj|, Semboloni 
et al. [35|), have demonstrated the potential of WL to 
place important constraints on the cosmological model. 
In particular, the CTIO results were among the first to 
indicate that the normalization of mass fluctuations was 
lower (cTg « 0.8) than the WMAPl (Spergel et al. (Hi) 
analysis suggested (as ~ 0.9), a finding confirmed now 
by the WMAP5 analysis (Komatsu et al. 

In this study, we concentrate on WL-detected clusters. 
WL identifies clusters as peaks in shear maps. The clus- 
ter abundance is one of the four most promising tools 
to measure dark energy (Albrecht et al. 01), together 
with supernova surveys, baryonic acoustic oscillations, 
and WL surveys. The cluster abundance is sensitive 
to the redshift dependence of the angular diameter dis- 



tance, the time dependence of the expansion rate and the 
growth rate of structure. There are three very appealing 
features of WL-detected clusters: 

1. The shear is straightforward to predict, given 
a model of structure formation. This is not 
the case for X-ray, Sunyaev-Zeldovich, or optical 
clusters, for which the conversion from measur- 
able quantities — flux, temperature, decrement — to 
virial masses depends on the complex gas physics 
of baryons in clusters. Cosmological constraints are 
likely to be biased if the mass-observable relation 
evolves with redshift in a way that mimics cosmo- 
logical changes. 

2. An immediate consequence of the first point is that 
shear peaks are easily studied from numerical sim- 
ulations of dark matter, since they are not very 
sensitive to the baryonic matter distribution. This 
is a huge advantage over the other cluster meth- 
ods, because dark matter simulations are relatively 
cheap to produce with the current computer tech- 
nology, whereas simulating the evolution of baryons 
depends on poorly known astrophysics and is com- 
putationally very demanding and expensive. 

3. Finally, shear peaks come "for free", i.e. a sur- 
vey designed to measure the "cosmic shear" 2-point 
function can detect clusters without much addi- 
tional observational effort. 

One issue that has raised concerns about the utility of 
WL clusters is the projection effect: the 2D shear maps 
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encode the 3D matter distribution information. If a clus- 
ter produces a detectable signal in a lensing map, then to 
this signal contributes not only the mass of the cluster, 
but also the matter along the line of sight. This mat- 
ter, external to the cluster itself, can be either nearby it, 
and therefore spatially correlated with it, or it can be at 
such a large separation that there is no correlation with 
the cluster. Hence, we can divide contributions from ex- 
ternal matter along the line of sight into correlated and 
uncorrelated projections. Thus, reconstructing the mass 
of a cluster can be challenging as there is no predefined 
way to separate its shear signal from that induced by 
other structures on the same line of sight. 

Many theoretical and numerical endeavors have been 
aimed at estimating the magnitude of this effect, its im- 
pact on cluster mass reconstructions, and on the cosmo- 
logical constraints derived from the WL-detected clus- 
ter abundance. Here we just mention a few of these 
works. Hoekstra [l^ used the aperture mass statistic, 
introduced by Schneider [3lj | , to evaluate the uncertainty 
in cluster mass estimates due to the uncorrelated large- 
scale structure (LSS) contamination. For large clusters 
at z = 0.1, they found the uncertainty to be « 6%. 

Other numerical studies have focused on the efficiency 
of WL surveys to detect clusters. Hamana et al. [ll| com- 
pared the signal-to- noise (S/N) of the peaks measured 
from simulated convergence maps to the Navarro- Frenk- 
White (NFW) (Navarro et al. |3) S/N of the 3D clus- 
ters producing the peaks. They found scatter and bias 
in the predicted vs. measured S/N relation, which they 
explained through both departures of the 3D halo den- 
sity profiles from the spherically symmetric NFW, as well 
as LSS projections. Tang and Fan [i^l also attributed 
a similar scatter to the asphericity of halos, and to halo 
substructure. De Putter and White [8(| concluded that in- 
dividual cluster shear measurements are contaminated by 
projections of very small clusters along the line of sight, 
i.e. halos with masses smaller than 10^^ Mq. Hen- 
nawi and Spergel performed an extensive numerical 
study of the WL-cluster detection. The authors used to- 
mographic information to determine the redshifts of the 
lenses, along with an optimal filter to enhance the cluster 
detection. One of their results was that even in the ab- 
sence of the intrinsic ellipticity noise - which is known to 
produce spurious peaks, see for instance Hamana et al. 

- about 15% of the most significant peaks were due 
to LSS projections, and not to any halo in particular. 
Dodelson and Maturi et al. [23] have independently 
proposed optimal linear estimators with a shape designed 
to filter out the LSS contributions. Maturi et al. [1^ used 
A'^-body simulations to compare their estimator, which 
is formally an aperture mass estimator, to the standard 
aperture mass filter, proposed by Schneider et al. [s^l 
for measurements of cosmic shear. The numerical tests 
have indicated their filter to be more efficient than the 
standard aperture mass at discarding 'false positives', i.e. 
those shear peaks arising from LSS projections. Schirmer 
et al. [lOl and Maturi et al. [23| detected clusters in shear 



data and found the performance of the optimal filter of 
Maturi et al. [1^ to be similar to that of the aperture 
mass statistic used with a NFW-matched filter. Another 
example of cluster detection in shear data is the work of 
Abate et al. [1]. 

Marian and Bernstein [2l| also used an optimal 
shear filter to forecast WL-cluster constraints for future 
ground- and space-based surveys. They found that cos- 
mological constraints are not severely altered by uncor- 
related LSS projections. One of the critical assumptions 
that forecasts adopt (see also Wang et al. ^l|) is to con- 
sider the Sheth-Tormen (ST) (Sheth and Tormen js^l) 
mass function a proxy for the shear-peak mass function. 
But is the shear-peak mass function as sensitive to cos- 
mology as the 3D one? In a pioneering study by Re- 
blinsky et al. [2§|, reasonable agreement was found be- 
tween the analytic predictions from the model of Kruse 
and Schneider [2C| and peak counts measured from nu- 
merical simulations. However, the small sample volumes 
studied, Vtot ~ lQ~^[h~^ Gpc]'^, limited the statistical 
power of the conclusions, especially for application to fu- 
ture WL data sets. In a recent analysis (Marian et al. 
[13]), we used a very large ensemble of numerical sim- 
ulations to pursue the issue further. We showed that 
the 2D mass function of peaks projected in thin slabs 
scales with cosmology in the same way as the measured 
3D mass function, and also as the ST prediction. Owing 
to the large total volume Vtot ~ l[/i~^Gpc]'^ per cos- 
mological model, the result was established at very high 
statistical significance. This suggests that for the pur- 
pose of constraining cosmology, one can avoid the difficult 
task of converting shear peaks into cluster masses. If the 
cosmology dependence of the shear peaks is understood, 
one needs only measure their abundance to obtain the 
desired constraints. Dietrich and Hartlap took this 
approach numerically and derived constraints from shear 
peaks measured from simulations with the aperture mass 
statistic. See also the related work of Kratochvil et al. 
[11], and Pires et al. [H]. 

In this paper we continue our investigation of the effect 
of correlated mass structures on the projected-peak mass 
function and its power to constrain cosmological models. 

To achieve this goal, we use numerical simulations. For 
each simulation box we divide the matter density field 
across the line of sight up into slabs. Our observable is 
the density field projected into such slabs. If the slabs 
are thin, this observable is equivalent to the WL conver- 
gence - defined as the surface mass density scaled by the 
critical surface mass density, Edit- Even when the slabs 
get thicker (256, 512 Mpc), the critical surface mass 
density varies only a few percents of its central value 
for typical lenses and sources in the ACDM model, as 
pointed out by Metzler et al. [2^ . The reason we are in- 
terested in slabs is that, by definition, we expect only the 
matter in a limited spatial volume around a main cluster 
lens to give rise to correlated projections. We underline 
that uncorrelated projections are the subject of a future 
study, and not of this study. 
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We extend the analysis of Marian et al. [24] in a num- 
ber of important ways. Firstly, we apply the halo model 
of structure formation to quantify the expected mass in- 
crease arising from correlated structures along the line of 
sight. Secondly, we generalize our earlier results through 
studying the projected peak mass function in slabs of 
various thicknesses. Thirdly, we show how the projected- 
peak mass function can be modeled. 

The paper is structured as follows. In section ijlll we 
present the numerical simulations that our study is based 
upon; in ijllll we use the halo model to predict the change 
in peak masses due to correlated projections and in PVI 
we describe the optimal filter that we apply to the data. 
In f|V]we show how we find and select the density peaks. 
Section Will contains the main results of our study, and 
in ijVIIII we draw our conclusions. Finally, we would like 
to clarify the mass definitions that we use throughout 
this work. We shall refer to those masses measured using 
full 3D information, as 3D masses (M3D). Similarly, 2D 
masses are measured from 2D data (M2d)- Depending on 
the context, M3D is usually either the virial mass defined 
by ST, or the mass defined by the friends-of- friends (FoF) 
algorithm of Davis et al. Q . We shall specify the meaning 
of the notation M3D in the context where we use it. M2D 
will always mean the mass of a projected peak, and the 
way we assign masses to projected peaks is described in 
section ^IVl 



II. NUMERICAL SIMULATIONS 



We have generated an ensemble of 32 A^-body sim- 
ulations for 4 different cosmological models, with 8 re- 
alizations per model. The simulations were run using 
Gadget-2 (Springel [1^), and they have 400^ particles 
of mass nip « 10^^ Mq in a simulated volume of 
512^ (/i-i M-pcf. The full details are in Table U There 
are 18 snapshots, logarithmically spaced in the scale fac- 
tor a, between redshifts 50 and 0. We used the 2LPT code 
of Scoccimarro [s^ to generate the initial conditions for 
the density fluctuations. For each of the 8 realizations, 
the initial conditions for the 4 cosmologies are the same. 
Thus, the cosmic variance on the comparison of measured 
quantities in different cosmologies will be minimal. The 
transfer function was taken from cmbfast (Seljak and 
Zaldarriaga [11] ) (the output at redshift 0). We identi- 
fied the simulation halos using the FoF algorithm, with 
b = 0.2, where b is the fraction of the inter-particle sepa- 
ration. The smallest of our halos have 30 particles. Our 
fiducial cosmology is described by Model 1 in Table [H 
We shall refer to the other models as the variational cos- 
mologies. 



TABLE 1: Numerical simulation details. Top-table columns 
are: matter, dark energy, and baryon density parameters; 
normalization of fluctuations; mass per particle in units of 
(10^^ Mq); number of realizations. Bottom-table columns 
show: box length; dark energy equation of state w — 
Pde/Pde; number of particles; Hubble's constant in units of 
100 kms~^ Mpc~^; baryon fraction. 
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Model 1 


0.27 


0.73 


0.0460 


0.9 


1.57 


8 


Model 2 


0.22 


0.78 


0.0375 


0.9 


1.28 


8 


Model 3 


0.27 


0.73 


0.0460 


0.75 


1.57 


8 


Model 4 


0.32 


0.68 


0.0545 


0.9 


1.86 


8 



L = 512 Mpc w^-l Np= 400^ h = 0.72 /b=17% 



III. THEORETICAL PREDICTION OF THE 
CORRELATED PROJECTION EFFECT 

In this section we develop a theoretical prediction for 
the correlated projection effect. We base our calculation 
for the correlated projections on the halo model, and we 
rigorously test it against numerical simulations. Our ap- 
proach is the following: we test the strength of the pro- 
jection effect in slabs of various thickness. For a cluster 
at the center of such a slab, we expect that the contribu- 
tion from structures correlated with it will flatten out as 
we increase the slab thickness. Only those objects within 
the correlation length of the cluster (ro ~ 30 Mpc 
for clusters more massive than 10"/i~1Mq) will affect 
its lensing signal. As the slab thickness becomes large 
- a few hundreds Mpc - it is reasonable to assume 
that chance or uncorrelated projections will start to con- 
tribute. 

If the large-scale distribution of the mass in the Uni- 
verse were perfectly described by the halo model, then 
all the mass would be contained in halos with universal 
properties: spherical symmetry, a universal density pro- 
file, and a mass function; see [5| for a halo model review. 
In this framework, we would like to estimate the average 
contribution to the projected density of a halo of mass 
Ml from structures correlated with it. The probability of 
finding a halo with mass in the interval (M2, M2 + dM2), 
within a range {r,r + d^r) of the center of Mi is given 
by: 

V{r,M2\Mi)dM2 d^r = ^(Ma) [l+^x {r, M2\Mi)]dM2 dV, 

a) 

where n{M2) is the halo mass function, i.e. the comoving 
number density of halos of mass M2 per unit mass, and 
is the cross-correlation between the center of halo Mi 
and the mass distribution of A/2. The expected density 
profile of matter external to the halo Mi is; 

(pext(r|Mi)) = J dM2 M2 n{M2) [1 + Cx (r, M2IM1)]. 

(2) 

The cross-correlation function can be obtained from the 
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FIG. 1: The change in mass due to correlated projections 
as a function of the Sheth-Tormen virial mass. The starred 
red symbols correspond to the optimal compensated filter de- 
scribed in section tjlVI and the empty green circles to the top- 
hat filter. These curves are for a slab thickness of 51 Mpc. 
The blue triangles show the maximum correlated projection 
effect for the compensated filter, as given by Eq. (|10p. This 
plot corresponds to the fiducial cosmology. 



cross-povirer spectrum: 



— fc2p^(fc,M2|Mi)jo(fcr), 



(3) 



where jo{x) = sin(a;)/x is the zero-order spherical Bessel 
function. The halo center cross-power spectrum, at lead- 
ing order in the halo perturbation theory, is: 



Px {k, M2IM1) = 6(Mi) b{M2) U{k, Ah)P{k) 



(4) 



where b{M) is the first order linear bias parameter, P{k) 
is the nonlinear matter power spectrum. J7(fc, M) is the 
Fourier transform of the 3D density profile of a halo with 
mass M: 



M 
MJo 



dr r'^p{r, M) ja{kr) 



(5) 



In the above, we have assumed that the 3D density pro- 
file of the halo is cut off at the virial radius i?vir- The ex- 
pected projected mass density external to the halo with 
mass Ml is the integral of Eq.(l2|) along the projection 
axis - for example the z direction: 



(Ee,t(rx|Afi)) =2 



L/2 



^i?.vir(Afi)2-r; 



_dz (Pcxt(2,r_L|Mi)) , 

(6) 



where r±_ is the 2D position vector transverse to the line 
of sight, and L is the slab thickness. For convenience, 
we have taken the center of halo Mi to be at half the 
slab's thickness. Combining the above expressions, we 
can write the last equation as a sum of contributions 
from uniformly-distributed matter, uncorrelated with the 
cluster Ml, and from clustered matter, correlated with 
the halo: 

(Mext(Mi)) = (A/uni(Mi)) -I- (Meorr(Afl)). (7) 

To predict the change in the mass of a projected peak due 
to both these terms, we need to convolve Eq. ^ with a 
filter, that we shall generically denote by W: 



(Afuni(Afi)) = 2p„, / d^r^W{r^) 



L/2 



dz : 



(8) 



(M,orr(Afi)) = 26(Mi) j dM2M2h{M2)n{M2) 

X j ^eU{k,M2)P{k) J d\^W{r±) 



L/2 



X / dz jo{k Jrl + z^) , (9) 

where pm is the mean matter density of the Universe. 

is a dimensionless function, so the right-hand side 
members of the above equations have mass dimension. 
If ly is a compensated filter (see section § IIV[) then the 
right-hand side of Eq. ([5]) is very small and negative; it 
would be in fact zero but for the dependence on rj_ of the 
lower limit of the line-of-sight integral. Throughout this 
analysis we are using a compensated filter, so we shall 
drop the contribution of (A/uni) since it is tiny compared 
to (Mcoii)- Note that an upper bound for the magnitude 
of the impact of correlated projections can be obtained 
if one takes the lower limit of the line-of-sight integral in 
Eq. ^ to be 0, and the slab thickness to be infinite. In 
that case, only the modes transverse to the line of sight 
contribute to (Mcon) and one can rewrite Eq. ^ in the 
following way: 

(M-,T(^A)> - b(Mi) I dM2M2b{M2)n{M2) 



X J dkj_k^U{k^,M2)P{k^) 

X / dr^r^W{r^)Mk^r^) , (10) 
Jo 

where Jq is the zero-order Bessel function. To illustrate 
the mass change induced by correlated projections on 
3D NEW halos, we have estimated Eq. Q in two cases: 
a top-hat filter, and a compensated, NFW-shaped filter 
that we describe in mV\ and that we used for the en- 
tire study. The limits of the /c-integral were adjusted to 
match the volume of our simulations. For the top- hat 
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filter, Eq. ^ can be written as: 



(Meo„(Mi)> = 2b{A'h)p^ / dkP(k) X 

Jo 

cos(fcz) - cos{k^/R~JM^Y+^) . (11) 



X / dz 

'flvir(Afl) 



Here we have taken U{k,M2) ~ 1; the integral over M2 
was then evaluated directly, giving p^. 

Figure [1] shows the change in mass arising from corre- 
lated projections in a 51 Mpc slab as a function of the 
virial mass of halos. The top-hat results are represented 
with green empty circles, while the red stars correspond 
to the compensated filter that we use throughout this 
work. The blue triangles are the upper limit contribution 
given by Eq. (jlOp . and computed with the compensated 
filter. The message of the figure is that the filter choice is 
crucial in reducing the impact of correlated projections 
on the measured peaks: the change in mass is « 8 times 
larger for the top-hat filter than for the compensated, 
NFW filter. This is due to the fact that the latter is 
very sensitive to cluster-scale modes, and less so to the 
larger-scale modes which give rise to projections. 



IV. FILTERING THE PROJECTED DENSITY 
FIELD 

In this section we present the filter with which we con- 
volve the density field, similar to that described in Mar- 
ian and Bernstein [l^l • We used this filter in section ^IIII 
to estimate the halo-model prediction of the mass change 
induced by correlated projections. Every point of the fil- 
tered density map has an associated mass given by: 



(12) 



where W is the filter function and E denotes the pro- 
jected density of dark matter. The filter is compensated 
[3lj, i.e. it vanishes when integrated over the aperture 
size. This ensures that no uniform mass sheet contributes 
to the measured projected density: 



d'^xWix) = 0, 



(13) 



where R is the aperture radius. The filter is also optimal, 
i.e. maximizes the signal-to-noise (S/N), and it is tuned 
to NFW clusters: an overdensity with the NFW density 
profile is detected with a maximum S/N. Although the 
profile of the halos best recovered by the filter is NFW, 
the mass is that defined by ST, i.e. an enclosed overden- 
sity threshold p = AvirPm as opposed to AvirPcrit- Here 
Pni is the mean matter density, pcrit is the critical density, 
and = f^mPcrit- Avir = 200 for ST and NFW. Thus, in 
the ACDM model, ST and NFW clusters have the same 
density profile, but ST halos have a larger cut-off radii 



and concentration parameters than NFW ones. The fil- 
ter is normalized to ST virial masses in order to facilitate 
the comparison of the 2D peaks with the 3D halos of our 
simulations, and also with the ST mass function. There 
is an agreement of « 10 — 20% between the FoF and the 
ST mass functions. As shown in numerous earlier studies 
[2ll . [23l . [bH , W is an optimal filter if it is proportional to 
the expected measured signal and inverse proportional to 
the variance of the noise. If the observable is the conver- 
gence field_(or in this case the projected density field), 
Schneider 31] has shown that W is optimal and compen- 
sated, if it has the form: 



W{x)=Cw (SsT(a;) -EsT(i?)) 



(14) 



with the mean density inside the aperture radius R de- 
fined by: S(i?) = 2/R'^ J^dxxT,{x). We determine the 
normalization constant Cw from the conditions that the 
filter be tuned to the NFW profile and return the ST 
virial mass: 



Jd^x\J:s^{x)\'~TTRl,^ tl^{R) 



(15) 



We choose the aperture radius R to be the virial ra- 
dius associated with an ST cluster of mass M^-J!, i.e. 
M^J = 47rAvir Pm ^st is the truncated ST pro- 

jected density: 

Est {x) = 2 6^'^ Pni /sT (x) . 

Vs is the scale radius of the cluster. The characteristic 
overdensity of the profile, S^"^ , is related to the concen- 
tration parameter Cst by the condition that the mean 
density within i?vir should be AvirPm: 



(5^^ — A • 



,/ (3 [log(l -f Cst) - Cst/(1 + Cst)]) ■ 



Finally, /st is a function which depends on cosmology 
only through the concentration parameter, see for in- 
stance Hamana et al. [11]. Its expression varies de- 
pending on whether it is evaluated at r < r^, r = rg, 
Tg < r < i?vir, and r > i?vir: 



2U/2 



cosh" 



(1-^^)(1 + Cst) 



-1 ( e 



+cst) 



1 



3(1 + Cst) 



1 



" Cst 
COS 



(l-e^)(H-CsT) 



-1 / _«!±£ST ^ 

V;(i+cst) J 



0. 



(16) 

Here ^ is dimensionless, ^ = r/rg. The concentration 
parameter is a function of the halo redshift and mass; 
we have followed the prescription of Bullock et al. Q to 
compute it and used the method of Smith and Watts 
[stI to map Cnfw Cst- The projected density field 
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FIG. 3: The unfiltered projected density of the same 51 h ^ Mpc slab shown in Figure (2] The scale of the picture is 130 x 
70/i~^Mpc. 
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FIG. 4: Same snapshot as above, but this time the projected density field is filtered. The detected peaks are marked with 
colored symbols depending on their mass range: purple squares for M > 5 x lO^** Mq; red circles for 5 x 10^* M© > 
M > 2 X 10" h-'^ Mq; black rombs for 2 x lO" Mq > M > 8 x lO" h'^ Mq; and orange ellipses for 8 x 10^^ h~'^ M© > 
M > 5 X 10^^ h-^ Mq. The FoF halos are represented with crosses and the color choice for each mass range is the same as for 
the peaks. The scale of the picture is 130 x 70 Mpc. 

maxima in the filtered map. A selection algorithm is then 
applied to decide on the validity and on the final mass of 
these peaks. 

One distinguishing feature of our method is that we 
filter the density maps recursively, with a hierarchy of 
matched filters. The shape of the matched filter depends 
upon the mass M^j'J of the cluster one is attempting to 
detect. We start with the highest-mass filter, convolve it 
with the map, and find all peaks above the target mass. 
We repeat the procedure for lower mass filters. Peaks 
that have been already found with a higher mass filter 
are discarded. Thus, we do not face the problem of 'halos- 
in- halos' (Hamana et al. flT'l). The last filter applied cor- 
responds to the detection threshold of 5 x 10^'^ h"^ M0. 
Below this threshold, the mass resolution of the simu- 
lations affects the filter efficiency, as shown in the next 
section. To the identified peaks we associate a virial ra- 
dius given by i?vir = (3M^;'^/47rAvirPm)^^'^, as described 
in the previous section. If two peaks have distinct centers 
(i.e. x — y coordinates), but overlapping radii of the same 
size, we keep both peaks. If the centers are distinct, but 
one of the peaks has a smaller radius than the other one, 
and it is therefore 'enclosed' in the projected disk of the 
larger peak, we keep the higher-mass peak and discard 
the lower-mass one. Figures [2J [3J and [5] show filtered 
and unfiltered maps of one of the slabs of the fiducial 



smoothed with the filter (|14p matched to a cluster of mass 
M^j^ will yield a maximum signal at the location of an 
ST cluster center, which we denote as M(a;o|M^;'^). If the 
cluster's mass were that of the filter, then the amplitude 
of the smoothed field at the peak would be: 

M{xa\Mll) = M^. (17) 

In practice, in order to assign a mass to a cluster, we 
apply several filters of different masses, and interpolate 
to find the mass M^j'J at which ((TT)) would be satisfied. 
This is the mass assigned to the peak. 



V. A HIERARCHICAL PEAK FINDING 
ALGORITHM 

As stated earlier, our observable is the projected den- 
sity field of dark matter, and our goal is to study the 
effect of correlated projections on the 2D-peak mass func- 
tion. We divide the simulation cube in several numbers 
of slabs of equal thickness: 20 (slabs of approximately 
26 ft-^Mpc), 10 [Ki 51/1-1 Mpc), 5 (w 102/1-1 Mpc), 2 
(= 256/1-1 Mpc), and 1 (= 512 /i"! Mpc). Each set of 
equal slabs is analyzed separately, i.e. the projected den- 
sity field is filtered, and peaks are identified as the highest 
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cosmology. Figure [2] presents a filtered slab, with the de- 
tected peaks and FoF halos marked on the filamentary 
structure of the projected density field. Figure [3] is a 
zoom-in of the same slab, before filtering, and Figure [4] 
is a zoom-in after filtering. In the latter, FoF halos are 
marked by crosses of different colors, according to the 
mass range where they belong, while peaks are marked 
by symbols of various shapes and colors - see the caption 
for the mass legend. 



A nontrivial issue regarding the slabs is that of the 
peaks which are situated right on the boundary of two 
slabs and are therefore split in two. The consequence of 
such a boundary split is that the peak is detected in nei- 
ther slab or, even if it is detected in one of them, its filter- 
associated mass is significantly reduced compared to the 
true value. On average, in a 50 h^^Mpc slab, about 
20% of the peaks suffer from boundary-split effects. We 
addressed this problem by considering two sets of slabs: 
'normal', and 'interlaced' slabs. Compared to the normal 
slabs, the interlaced ones are shifted by half a slab thick- 
ness on the axis of projection. Thus, the boundary of two 
normal adjacent slabs goes right through the middle of 
an interlaced slab, so that halos that are divided in the 
normal slabs appear whole in the interlaced slabs, and 
vice versa. We analyze the peaks in both sets of slabs 
and compare them based on their location in the x — y 
plane (we consider the projection axis as the z direction), 
and their filter-associated mass. To be precise, the peaks 
in every normal slab are compared to those in the two 
adjacent interlaced slabs which contain the normal slab. 
The peaks in every interlaced slab are also compared to 
those in the two adjacent normal slabs which contain the 
interlaced slab. If we find peaks with the same or nearly 
the same x — y coordinates, but with masses that differ 
by more than 5%, we consider that we have a boundary 
split case and select the peak that has the higher mass. 
If the difference in mass is smaller than 5%, we keep the 
original masses. 



Halos often 'live' in filaments; if the filaments are ori- 
ented along the projection axis, they can contribute to 
the mass of peaks sourced by the respective halos, or they 
can even form peaks of their own. Even when a halo is 
not right on a slab boundary, but it is close to it, its 
filament can extend beyond the boundary. By splitting 
the filament, the projected mass of the halo is slightly 
decreased. The 5% threshold mentioned above not only 
removes the problem of split halos, but it also makes split 
filaments less likely. For the thickest slab with potential 
boundary problems, i.e. 256 Mpc, we use a higher 
threshold of 15%, in order to avoid biasing our mass mea- 
surements to higher masses. In this case there is only one 
boundary, and hence fewer split objects. 



VI. TESTING THE FILTER 
A. Numerical tests of the filter 

The filter described earlier is optimal in the ideal con- 
ditions of infinite particle-mass resolution. The fact that 
the simulation halos have a discrete mass distribution 
changes the performance of any filter when applied to 
simulation data. Also, our filter is tuned to spherically 
symmetric halos, but real halos are more triaxial than 
spherical. We have investigated the efhciency of our fil- 
ter in two idealized cases: first, on synthetic NFW halos, 
to understand the effects of the finite mass and mesh 
resolutions; second, on isolated FoF halos (that is, only 
the particles from a single halo) of our simulations, to see 
the effects of the halo asphericity. Both tests were carried 
out for the fiducial cosmological model, as there is no rea- 
son to assume that such numerical effects would manifest 
themselves differently in the variational cosmologies. 

The synthetic NFW halos that we have generated have 
the same mass resolution as our simulations. The par- 
ticles are sampled in the radial direction with the NFW 
density profile, and their angular distribution is uniform 
random on the sphere. The tests were performed for 3 
different masses, using an ensemble of 1000 synthetic ha- 
los for each mass. We found that the finite particle mass 
resolution mostly affects the signal in the inner region 
of the clusters, defined by the scale radius. That is be- 
cause the NFW density profile has a significant fraction 
of its mass enclosed within the scale radius (« 20%). 
Given the mass resolution of our simulations, the inner 
region of average clusters is poorly populated with par- 
ticles, and actually unresolved for radii smaller than the 
softening length of the simulations (60 kpc, in our 
case). Therefore, we have chosen to 'core' the ST pro- 
file, i.e. we consider the density profile of the filter to be 
constant for a radius < rcorc^ 

V (^\ — J ^ST(''corc) , f l£ ''core f^Q\ 
I ^stI^' ) 1 1^1 core 

Thus, the traditional central cusp of the ST profile has 
been removed from the filter, and the signal of the filtered 
halos is redistributed so that the outer layers contribute 
more to it. Even if our simulations had allowed us to 
use the ST profile without removing its cusp, there is 
a more physical reason that justifies our choice. Halos 
relevant to WL surveys have an Einstein radius which is 
only 2-3 times smaller than the scale radius. Since we 
are interested in the WL regime, the measured signals of 
the halos must come from regions well outside the Ein- 
stein radius, or else the WL approximation breaks down. 
Note that the filtered density will be estimated from real 
data by applying a corresponding filter to the shear map. 
The shear filter corresponding to our constant-core ST 
profile has a 'hole', i.e. discards the signal coming from 
within the core radius. This is necessary in practice since 
it is very difficult to measure the shear around clusters 
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FIG. 5: The dependence of the filter efficiency on the 
core radius, for the Fourier mesh used throughout the pa- 
per (60 kpc). The efficiency is defined as M2D/iVf3D, 
where AfsD is the Sheth-Tormen virial mass of the syn- 
thetic halo, and M2D is the mass of the projected peak, 
assigned as described in ijlVI We show the average of 
measurements of 1000 projected synthetic Sheth-Tormen 
halos, in the fiducial cosmology. We chose a core radius 
of 200 kpc to perform our analysis. 
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FIG. 6: The dependence of the filter efficiency on the 
Fourier mesh, for a core radius of 200 ft"^ kpc. The fifier 
efficiency is defined as in Figure [S] For each mass we plot 
the average of the measurements of 1000 synthetic Sheth- 
Tormen halos in the fiducial cosmology. Throughout our 
analysis, we used a mesh of 60 kpc. 



at small radii from the center: most shear measurement 
methodologies fail in regions of strong lensing, and clus- 
ter member galaxies often obscure the lensed background 
galaxies. 

Figure O plots the dependence of the filter efficiency 
- defined as Af2DA^3D - on the filter core radius for an 
ensemble of synthetic halos of mass 10"'^"' h^^ M©. M2D is 
the mass that we assign to the density peaks in the way 
described in section ijlVl and M^d is simply the ST virial 
mass of the synthetic halos. The figure clearly shows 
the filter efficiency to be « 90% even for halos obeying 
perfectly the NFW profile, if no coring is done. Though 
not shown in the figure, we have also probed that the 
efficiency would be more than 95% if the mass of the 
simulation particles were lower by a factor of 10, with no 
core applied. However, since this mass resolution is not 
available, throughout this study we have adopted a core 
value of 200 h^^kpc. This is also valid for the results 
published in Marian et al. |22| . 

Equation p2| is a real space convolution. By smooth- 
ing the density field of dark matter with a filter, we ex- 
pect to erase the structure on scales smaller than the 
scale of the filter. It is convenient to do the smoothing 
in Fourier space, and in order to compute the discrete 
Fourier transform of the density field, we project it on a 
grid. We refer to the ratio of the simulation cube side 
L and the number of grid points taken on that side, as 



the Fourier mesh resolution. Our whole analysis of the 
simulation data was done with a mesh resolution of 60 
/i~^kpc. In Figure [6] we plot the filter efficiency depen- 
dence on the Fourier mesh resolution for the chosen core 
value. Large halos are less affected than small halos by 
resolution effects, and in general the filter efficiency is of 
about 98% for our mesh and core chosen values. Since 
we do not wish to include in our analysis halos that are 
too small and too 'poor' in particles, we have imposed a 
mass detection threshold of 5 x 10^"^ h^^ Mq for the pro- 
jected peaks. We do not include in our analysis peaks 
with mass below this threshold. 

Figure [7] shows the correspondence of FoF masses and 
2D masses recovered by our filter, in the absence of pro- 
jections. The FoF halos have been projected by them- 
selves: only the particles identified by the halo finder as 
being part of a halo have been projected. Again, the effi- 
ciency is defined as M2D /Msu : M2D is the mass assigned 
to the density peaks according to section ^IVl and A/sd 
is the FoF mass of the simulation halos that we analyze. 
The filter efficiency is depicted for 3 large mass bins. The 
FoF halos were selected from all eight realizations of the 
fiducial cosmology. The plot reveals that the vast major- 
ity of the FoF halos are shifted to lower masses by our 
filter. This is mostly due to departures of the halo den- 
sity profiles from the spherically symmetric NFW profile: 
the filter given by Eq. is a spherical overdensity (SO) 
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FIG. 7: The filter efficiency in tfie absence of projections. AfsD is tlie FoF mass of tlie simulation halos, and Af2D is the mass 
of the peaks obtained by projecting the FoF halos by themselves, i.e. only the particles forming them. Three FoF mass bins 
are considered, with the total number of halos per bin indicated in the figure. We attribute the shift towards lower masses to 
departures of the FoF halos from the spherically-symmetric NFW profile. The dispersion in the recovered mass decreases for 
the bins of lower mass, in accord with the fact that the largest halos are the most triaxial ones. 



filter that we apply to FoF halos. FoF and SO masses 
are known to differ. Indeed, had we applied an SO 3D 
halo finder, or had we used another FoF linking length, 
the filtered masses would have had a different distribu- 
tion and mapping to the 3D halos. It is also possible 
that the scatter in the concentration parameter, which is 
well known to be described by a log-normal distribution 
(Bullock et al. fsj) might also cause some of the scat- 
ter seen in Figure [71 Also, halos have substructures, e.g. 
[1^. Since we are interested in the cosmology scaling of 
the projected-peak mass function, and not in the analy- 
sis of individual clusters, the mapping between FoF halo 
masses and filtered masses is not, in this study, of great 
significance to us. After all, there is no absolute way to 
define a halo mass. It is worth emphasizing that, unlike 
Metzler et al. [1^ and Clowe et al. Q our recovered 2D 
masses are not the result of averaging projections of the 
same clusters along different lines of sight. The halos are 
projected on the same axis as the simulation cubes, and 
to each of them we apply our filter 'blindly', i.e. without 
assuming any prior knowledge of the true position of the 



cluster center. This is valid for our whole analysis. Had 
our individual cluster mass estimates been an average 
over several lines of sight, we expect to have found a re- 
duced scatter and a smaller skewness in the distributions 
shown in Figure [71 

Figure [51 shows the correspondence of FoF masses and 
2D masses in the presence of projections. Here we have 
matched the FoF halos of the simulations with the peaks 
projected in slabs of various thickness. The different col- 
ors indicate the thickness of the projection slabs. The 
binning of the FoF masses is finer than in the previ- 
ous figure. The sharp cut-off visible in the right low- 
est panels stems from our selection of peaks larger than 
5 X 10^^ Mq. The halos sourcing these peaks can 
however have masses lower than this threshold. Not sur- 
prisingly, the same trend present in Figure [71 is also visi- 
ble here. Again the FoF halos are shifted towards lower 
masses, and the distribution of peaks is non-gaussian. 
A log-normal distribution fits some of the intermediate 
mass bins. The slab thickness does not seem to play 
a significant role in shaping the distribution. While in 
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FIG. 8: The filter efficiency in tfie presence of projections. Tlie efficiency is defined as M2u/Msu, wliere M^u is the mass of 
the FoF halos that source the peaks measured in different slabs. A-hu is the mass of the peaks formed by the FoF halos, when 
projected together with the rest of the matter in the slabs: 26 Mpc - magenta- , 51 Mpc -cyan- , 102 Mpc - red 
- , 256 h-^ Mpc - green - , 512 h ^ Mpc - blue. The distribution is non-gaussian, and the spread is mostly due to departures 
of the 3D halos from the spherically-symmetric NFW profile. 



this study Figure [5] is just illustrative, should we wish 
to predict the projected-peak mass function from the 3D 
one, we would need the distribution M2D ~ ^^3D, as sug- 
gested by Metzler et al. [1^. Figures [7] and [5] serve as a 
warning against blind comparisons between 3D and 2D 
quantities, as often seen in the literature. Figure [7] is an 
example of 'filter effects', i.e. scatter between 2D and 3D 
masses due to different mass definitions and other causes. 



but completely unrelated to LSS projections. Figure [5] 
shows that such filter effects are dominant compared to 
the correlated projection effect. Unless the 2D and 3D 
peaks are found with the same filter, it is very difficult 
to measure the impact of correlated projections directly 
on 3D masses. In order to remove the filter effects shown 
in Figure [71 one should measure the change in the mass 
of peaks projected in slabs of different thickness. This is 
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FIG. 9: The fractional difference of the number of projected peaks in slabs of 25, 102, 256, 512 Mpc and the peaks in the 
fiducial slab of 51 h^^ Mpc. Each panel corresponds to one cosmology. The difference in the mass functions for different slabs 
is due to the correlated projected term described by Eq. Q, as well as to uncorrelated projections in the case of the thickest 
slabs. 



how we proceed in the following sections. 



VII. RESULTS 

In our previous study Marian et al. (2^ . we ana- 
lyzed the effect of correlated projections in slabs of 51 
Mpc. Here we extend our analysis to several slab 
thicknesses: 26 h'^ Mpc, 102 h'^ Mpc, 256 h"^ Mpc, and 
512 /i-^Mpc. 



A. Slab thickness tests 

First we probe the dependence of the projected density 
peak abundance on the slab thickness. Figure [9] shows 
the fractional difference of the peak number counts in 
the slabs that we investigate, compared to our 'fiducial' 
slab of 51 Mpc. We notice the fractional differences 



deviate from 0, particularly at the high-mass end, an an- 
ticipated effect of projections. The deviations depend on 
the slab thickness and also on the cosmological model 
- the low-cTg cosmology shows greater sensitivity to the 
slab variations than the other three models. We expect 
the most significant changes from the fiducial slab peak 
abundance to occur precisely at the high-mass end: there 
are few high-mass clusters, so a small variation in their 
mass can alter their distribution on the exponential tail 
of the mass function considerably. In particular, for the 
512 Mpc and even the 256 Mpc slabs, this varia- 
tion may be increased by the presence of chance projec- 
tions along the line of sight, which seems to be suggested 
by Figure [SI The correlation length of large clusters is 
w 30 Mpc, therefore we do not expect correlated pro- 
jections to affect the abundance of peaks in these slabs. 
There is also an increase in the high-mass peak abun- 
dance corresponding to the fiducial slab compared to the 
26 /i~^Mpc one, while only a very small rise is visible 
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FIG. 10: The fractional change in mass due to correlated projections for peaks measured in different slabs and the fiducial 
slab peaks. Every peak projected in slabs with the fiducial thickness, was then tracked among the peaks projected in the other 
slabs. 



in the 102 Mpc slab compared to the fiducial one — a 
sign that the correlation function of clusters drops signif- 
icantly at scales larger than 50 Mpc. Overall, we are 
satisfied with the robustness of the mass function pro- 
jected in slabs of various thickness: analyzing the pro- 
jected density field with a compensated filter has ren- 
dered the resulting peak abundance largely independent 
of the slab thickness. 



B. The mass change induced by correlated 
projections on individual peaks 

In section ijllll we used the halo model to predict the 
change in the projected mass of a cluster induced by 
structures correlated with it, i.e. Eq. We have eval- 
uated this expression for the particular case of our filter, 
described by Eq. (|14p. In this case, Figure [1] indicated 
that there should be a very small change (< 2%) in the 
2D mass of peaks - defined in section ijlVI - , compared 



to their 3D mass - defined as the ST virial mass. In this 
section we test the halo model prediction. Since our 3D 
masses have not been identified with an SO filter, but 
with the FoF algorithm, we cannot probe the result of 
Eq. ^ directly. Instead, we circumvent this difficulty 
by comparing only masses of objects identified with the 
same filter. 

First we follow the mass change for the simulation clus- 
ters when projected in slabs of various thickness. We 
have traced every peak in the fiducial slab among the 
peaks corresponding to the other slabs. We have selected 
only those peaks that maintained roughly the same x — y 
coordinates in all the slabs where they were projected - 
we did not allow the coordinates to shift by more than 
« 120/i~^kpc. This restriction was imposed to ensure 
that we compare the mass of the same peaks traced in 
different slabs. Figure [10] presents the fractional change 
in mass for the peaks projected in the variational slabs 
compared to the fiducial slab (symbols). The error bars 
represent errors on the mean of the eight realizations of 
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on our data for a direct comparison with the halo model 
prediction shown in Figure [TJ The result is shown in 
Figure 1111 Some caution is necessary when interpreting 
this figure, as for some of the intermediate mass bins, the 
ratio M2T>/M^Y) goes below unity. Naturally, this should 
not be the case. We explain this anomaly by the fact that 
when an FoF halo is projected together with the struc- 
ture in which it is embedded, the surrounding matter 
renders it more homogeneous and the compensated filter 
may therefore attribute it a smaller mass than when the 
halo is projected by itself. However, despite this anomaly, 
FigurefTTlshould be taken as supportive of the halo model 
prediction, i.e. Eq. ([9]). Hence we expect that correlated 
projections alter the masses of clusters at the percent 
level, and therefore are not a dominant systematic for 
WL cluster mass measurements. If there is a projection 
bias in WL mass determinations, it is more likely due to 
structures uncorrelated with the main lens. 



FIG. 11: A test of Eq. M30 is the mass of FoF halos 
projected by themselves, and filtered according to Eq. (|14|l . 
M2D is the mass of the peaks sourced by the same FoF halos 
when projected in slabs of the fiducial thickness, together with 
all the mass in the respective slabs. The solid line is the halo 
model prediction from Eq. (|9}. 

each cosmology. We have measurements for all four cos- 
mologies. One remarkable feature may be noted: the 
simulation data indicate the impact of correlated projec- 
tions on the measured halo masses to be very small, i.e. 
about 2% for peaks measured in slabs of 26 Mpc and 
102 Mpc relative to the fiducial slab peaks. The ex- 
ception is the low-ijs cosmology, which shows stronger 
deviations at the high-mass end. The halo model predic- 
tion for the change in mass of peaks projected in slabs 
of 26 Mpc and 102 Mpc compared to peaks pro- 
jected in slabs of 51 Mpc is that it should be smaller 
than 0.1 %. FigurefTOlindicates a more noticeable change, 
however. The slightly larger measured mass changes 
( 4%) for the two thickest slabs are likely due to chance 
projections. 

Second, we carry out a comparison between peaks mea- 
sured in slabs of thickness and peaks measured in the 
fiducial slab. To be precise, we take a large number of 
the FoF halos 8000) from the fiducial cosmology, and 
project only the particles within the halos themselves. 
An identical procedure was followed to obtain the dis- 
tribution in Figure [T] In this case, we filter a larger 
number of halos, and take finer mass bins. We consider 
the mass of the peaks obtained through the projection 
of the FoF halos by themselves as a proxy for the 3D 
ST virial mass. We then compare this mass to the mass 
of the peaks sourced by the same FoF halos when pro- 
jected together with all the matter in slabs of the fiducial 
thickness. We do not expect this approximation to be 
perfect, but it is the closest measurement we can make 



C. Predicting the projected mass function 

We now try to address two of the main objectives of 
this work: in the context of correlated projections, we 
examine the behaviour of the 2D mass function and also 
the possibility to predict this behaviour. In Marian et al. 
[2^ we have shown that the scaling with cosmology of 
the 2D mass function projected in slabs of 51 Mpc 
follows very closely the scaling of the 3D mass function 
and also that of the semi-analytical ST model. 

Figure [T^ is similar to Figure 2 from Marian et al. [23 |. 
but generalizes that result to other slab thicknesses. The 
figure shows the fractional difference of the peak counts 
for the variational cosmologies scaled by the fiducial cos- 
mology peak counts (solid points). It also shows the same 
fractional differences for the 3D FoF mass functions mea- 
sured from the simulations (open magenta symbols), as 
well as the ST predictions (solid lines). The error bars 
are on the mean of the eight realizations of each cos- 
mology, and the cosmic variance between realizations is 
minimized by the choice of initial conditions mentioned in 
section f|TTl The difference in the SO filter and FoF mass 
definitions is also minimized: we compare 2D peaks with 
2D peaks and FoF halos with FoF halos. 

We maintain the conclusion of Marian et al. re- 
garding the very similar scaling with cosmology of the 
2D, 3D, and ST mass functions. 

We next try to fit the lensing mass function corre- 
sponding to the slab thickness of 26 Mpc. We shall 
assume that the functional form proposed by [s^ can be 
used to describe the projected mass function, but with 
the parameters of the model recalibrated. The ST mass 
function is given by: 

dnsT{M) Pm , , , . 

^l-^rflogM=^/s.(.)d., (19) 
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FIG. 12: The mass function of the variational cosmologies, scaled by fiducial cosmology mass function, considering slabs of 
different thickness. Filled symbols are peak measurements, open ones are FoF halos, and the continuous lines indicate the 
Sheth-Tormen mass function. The error bars represent errors on the mean of 8 realizations of each cosmology. 



with v = Sl^jcP' and: 

h-.{v) = [1 + (« [-civl2) . (20) 

In the above dn^-r = ?t.st ij^l^ AM gives the number den- 
sity of halos, is the variance of the top-hat-smoothed 
linear density field, 5sc = 1-69 is the critical density for 
spherical collapse. The ST model has two free param- 
eters, with the following values: a = 0.707, p — 0.3. 
The normalization parameter A is obtained from the con- 
straint that all mass in the Universe is in halos: 



extended Poisson likelihood for our simulations: 



dMMnsT{M) = p. 



(21) 



This leads to: A{p) = [1 + r(l/2 - p)] . For 

the ST mass function, A = 0.3222. The projected mass 
function fit conserves this normalization relation, but 
changes the values of a, p. To fit t he p arameters we fol- 
low the approach of Warren et al. [44I and compute the 



InA(q) 



EE 

1 = 1 i=l 



In- 



(22) 

where N-^ is the number of realizations per cosmology - 
8 in our case - , and N,^ is the number of mass bins con- 
sidered, riij is the measured number of projected peaks 
in mass bin i of realization j , while /i^ is the predicted 
number of halos for mass bin i and realization j. We 
find the point in the parameter space defined by q—{a,p) 
which maximizes the above likelihood, and thus we de- 
termine the best fit for the projected mass function. Note 
that, had we been attempting to constrain the mass func- 
tion at lower masses, we expect that the Poisson model 
would have not captured the true variance in the dis- 
tribution, owing to sample variance fluctuations (Crocce 
et al. Hu and Kravtsov [IBl)- 

We have fitted the free parameters for each cosmolog- 
ical model, and our results are the following: 
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FIG. 13: Fractional difference of measured and fitted 2D 
mass functions for all four cosmologies. The slab thick- 
ness is 26 Mpc. 



FIG. 14: Fractional difference of the 2D mass function 
measured in slabs of 26 Mpc and the prediction given 
by Eq 



1. Model 1: a = 0.719, p = 0.298, A = 0.3242 

2. Model 2: a = 0.690, p = 0.306, A = 0.3161 

3. Model 3: a = 0.665, p = 0.318, A = 0.3034 

4. Model 4: a = 0.727, p = 0.292, A = 0.3301 

Figure [T51 shows a comparison between the measured and 
the fitted projected mass functions for all cosmologies. 
We note that the ST functional form is not ideal to fit 
the high-mass end of the 2D peak function: the accuracy 
of the fits is better than 10% for the low and intermedi- 
ate mass bins, but drops to 30%-40% for the high mass 
bins. This is just a consequence of having projection 
noise added to measurements on a steep mass function. 
We also try to predict the projected mass functions of 
the variational models by using only the fit to the fiducial 
cosmology and ratios of the standard ST mass functions: 



fiducial 
"fit ^ 



n. 



X 

ST 



T fiducial * 



(23) 



In the above, n-^ is the predicted 2D mass function for 
any of the variational cosmologies, n^'^^'^^'^^ is the fit for 
the 2D mass function of the fiducial cosmology; and 
nf^""'^' are the standard ST mass functions for the vari- 
ational and fiducial cosmologies respectively. 

Figure [HI presents the comparison between this predic- 
tion and the measured 2D mass functions. The low- and 
high-f^m cosmology mass functions are captured with an 
accuracy of 10%-20%, at low and intermediate mass bins, 
while the Iow-cts model is worse. The modeling fails com- 
pletely at the high-mass end, which was to be expected 



given that for this mass range: (i) there is a large inac- 
curacy in the fit of the fiducial cosmology; (ii) the ST 
functional form does not describe well the 2D mass func- 
tions of the variational cosmologies either. This simplis- 
tic modeling should just be taken as a proof of concept, 
that the scaling with cosmology of the projected mass 
function presented in Figure [T^] can be used as a start- 
ing point for more complicated models of the 2D peak 
abundance. 



VIII. CONCLUSIONS 

Future or near-term WL surveys will be able to detect 
clusters as peaks in shear maps. Projections from LSS 
are likely to make it difficult to convert the amplitude of 
the shear peaks to traditionally-defined cluster masses. 
However, it is possible to circumvent this conversion if 
the cosmology dependence of the shear-peak counts is 
understood. In this case, one can constrain cosmology 
by measuring directly the shear-peak abundance. While 
this can be done in dark matter simulations, as recently 
shown by Dietrich and Hartlap [3], it is not practical 
to rely only on numerical methods. The standard cos- 
mological model contains 6 free parameters, and so the 
space of models to be tested may quickly become large. 
Thus, if a well-understood analytical framework can be 
identified for the peak counts, it will hopefully facilitate 
an accurate interpolation between expensive numerical 
predictions made in this high-dimensional space of cos- 
mological models. This is the main motivation for our 
study. 
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In this paper we focused attention on estimating and 
understanding the effect of 'correlated' projections on the 
mass function of density peaks projected in slabs of differ- 
ent thicknesses. We reserve 'uncorrelated' projections for 
future study. Correlated projections are ideally studied 
in slabs, given the fact that the correlation length of aver- 
age clusters is ~ 30 Mpc. In such slabs, the projected 
density is equal to the convergence up to a constant, so 
it is perfectly adequate to choose it as an observable. 

In f|ll] we described the large ensemble of numerical 
simulations that were generated for our study. These 
consisted of 4 cosmological models with 8 simulations 
per model and with total volume per model Vtot = 
l[/i-iGpc]3. 

In mil] we presented a calculation based on the halo 
model of structure formation of the expected change in 
the mass of 3D clusters induced by correlated projections. 
We showed that the mass increase depends on the filter 
applied to the projected density field. For a top-hat filter, 
the increase in mass due to correlated projections is of 
10 — 15%. For a matched, optimal, compensated filter, 
it is roughly a factor of 8 lower. Thus, the halo model 
predicts that correlated projections should not alter the 
virial mass of clusters by more than 2%. 

In [|V] we described our algorithm for identifying 
projected-density peaks. One important new element 
was that we offered a solution to the problem of 'peaks- 
in-peaks'. We did this through applying a hierarchical fil- 
tering strategy. We used different scales for the matched 
filter and took a top-down approach to peak detection. 
Smaller peaks that were contained within higher mass 
peaks were discarded. 

We have used the same optimal filter as in Marian and 
Bernstein |2l]|. This filter is compensated, maximizes the 
S/N, and matches the NFW density profile. Its normal- 
ization is chosen so that it returns the ST virial mass 
of halos. Tests on synthetic data, i.e. perfect NFW ha- 
los, have shown in section WI Al that the filter requires 
'coring' when applied to simulated data, due to the finite 
mass and mesh resolutions. The coring will also be nec- 
essary for application to real WL measurements, due to 
practical difficulties of making precise WL measurements 
in the inner regions of galaxy clusters. 

We have measured the distribution M2D /M^d - where 
M3D = M-poF ~ in the presence of projections: we have 
matched the FoF halos of the simulations and the peaks 
measured from the density field projected in slabs of dif- 
ferent thickness. We have found this distribution scat- 
tered and biased towards values smaller than 1. We have 
also measured the same distribution in the absence of 
projections, by projecting the FoF halos by themselves. 
Even in this case the distribution M2D/Af3D was scat- 
tered and biased towards values smaller than 1. This 
shows that most of the scatter between M2D and M3D is 
not due to LSS projections, but rather to the departures 
of the FoF halos from the spherically-symmetric NFW 
profile, to the existence of substructures, and possibly to 
the scatter and stochasticity of the concentration param- 



eter. 

The point of these tests is to clarify a somewhat loose 
perception in the literature concerning the impact of LSS 
projections on cluster mass measurements: when com- 
paring 2D versus 3D quantities (i.e. mass, S/N, etc.), 
one must account for the halo/cluster identification al- 
gorithms. The scatter in the relation between SO and 
FoF-identified halos is sure to confuse the interpretation 
of LSS projections, and can lead to an overestimation of 
the latter. Since we were unable to apply our matched 
filter algorithm in 3D (we used a 8192^ FFT in 2D, and 
so would require a 8192'^ FFT in 3D), we avoided the 
confusion of SO-FoF by examining the change induced 
by correlated projections in density peaks projected in 
slabs of various thicknesses. 

In Will we measured the projected-peak function in 
slabs of 26, 51, 102, 256, 512 h-^Mpc. The use of the 
compensated filter rendered the 2D mass function quite 
stable with the slab thickness. We then measured the 
mass evolution of peaks when projected in the different 
slabs. We found that the average change in mass of peaks 
measured in slabs of 26 h^^ Mpc and 102 Mpc, com- 
pared to their mass in the 51 Mpc slab, is no greater 
than 2%. For the thicker slabs of 256 h'^ Mpc and 512 
h^^ Mpc, the average change is also less than 4%; we 
expect these thicker slabs to be affected by uncorrelated 
projections. The halo model suggests that a peak mea- 
sured in a 51 Mpc slab should suffer a change in mass 
smaller than 0.1% compared to when it is measured in 
a 26 Mpc slab. So these changes are larger than 
the halo model predictions. On the other hand, we have 
tested that the absolute change in the mass of halos is also 
at the percent level, in accord with the halo model results. 
Let us emphasize that all these findings depend on the 
choice of filter: our filter reduces dramatically the impact 
of correlated projections compared to other filters. The 
importance of the filter function in reducing the LSS con- 
tamination of WL cluster mass measurements was first 
suggested by Maturi et al. [1^. We conclude that with 
a proper choice of filter, correlated projections are not a 
major contaminant for WL-cluster mass measurements. 
Uncorrelated projections are likely to be the dominant 
contributor to the projection bias. 

We also confirm the main result of our earlier study 
Marian et al. [2^: the projected peak abundance has 
the same scaling with cosmology as the 3D and ST mass 
functions. This is valid for all the slab thicknesses that we 
probed. We proposed a very simple fitting procedure for 
the 2D mass function, adopting the functional form of the 
ST prediction. Following the work of Warren et al. [43 |. 
we have estimated the free parameters of the ST model, 
by maximizing the extended Poisson likelihood. For all 
cosmologies, the fits have an accuracy of up to 10%-20% 
at low and intermediate mass bins. We also used only 
the fit for the fiducial cosmology 2D mass function, to 
predict the 2D mass functions for the other cosmologies. 
In this case, the accuracy was of about 20% in the same 
mass range. 
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Our goal for the future is to extend our analysis to 
uncorrelated projections, and to study both thcjoretical 
and numerical avenues for the prediction of the shear- 
peak abundance, and its cosmology dependence. Should 
such an undertaking be successful, it will establish the 
WL method for cluster detection as a very reliable way to 
extract cosmological constraints, limited only by survey 
systematics. 
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